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Abstract 

The following effects in the nearly forward ( "soft" ) region of the LHC are proposed 
to be investigated: 

• At small |t| the fine structure of the cone (Pomeron) should be scrutinized: a) a 
break of the cone near -0.1 GeV^, due to the two -pion threshold, and required 
by t-channel unitarity, and b) possible small-period oscillations between t = and the 
dip region. 

• In measuring the elastic pp scattering and total pp cross section at the LHC, 
the experimentalists are urged to treat the total cross section at, the ratio p of real 
to imaginary part of the forward scattering amplitude, the forward slope B and the 
luminosity C as free parameters, and to publish model-independent results on dN/dt. 

• Of extreme interest are the details of the expected diffraction minimum in the 
differential cross section. Its position, expected in the interval 0.4 < — t < 1 GeV^ at the 
level of about 10~^ mb- GeV~^H-10~^ mb- GeV~^, cannot be predicted unambiguously, 
and its depth, i.e. the ratio of da/dt at the minimum to that at the subsequent 
maximum (about —t = 5 GeV^, which is about 5) is of great importance. 

• The expected slow-down with increasing \t\ of the shrinkage of the second cone 
(beyond the dip-bump), together with the transition from an exponential to a power 
decrease in —t, will be indicative of the transition from "soft" to "hard" physics. 
Explicit models are proposed to help in quantifying this transition. 

• In a number of papers a limiting behavior, or saturation of the black disc limit 
(BDL), was predicted. This controversial phenomenon shows that the BDL may not 
be the ultimate limit, instead a transition from shadow to antishadow scattering may 
by typical of the LHC energy scale. 
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1 Introduction 



In this paper, some crucial issues that may be useful in preparing the experiments at the 
LHC are discussed and clarified. 

Measurement of the total proton-proton cross section will be one of the first priorities at 
the LHC. The importance of these measurements is two-fold. First, they are mandatory to 
fix the normalization of all subsequent measurements. Furthermore, the value of the total 
cross section will drastically narrow the range of the existing models with predictions for 
the total cross section ranging as crtot(14 TeV) = (125 ± 25) mb [1,2] or even more^. The 
knowledge of the total cross section will help in selecting a class of models of diffraction, 
based on the dominance of multi-gluon, or Pomeron (IP) exchange. 

Diffractive events, for example, diffractive Higgs production, are widely believed to pro- 
duce the cleanest signal of possible new phenomena [3-7]. 

Quantum ChromoDynamics (QCD), complemented with the Regge pole theory, and the 
unitarity condition superimposed, form the theoretical basis of the strong interaction. Both 
QCD and the Regge pole theory need experimental verification to clarify the role of higher 
QCD corrections, on the one hand, and to restrict the existing fiexibility in the Regge pole 
models, on the other hand. The scattering amplitude at the LHC can safely be parameterized 
by the dominating Pomeron exchange, appended by a possible tiny Odderon contribution, 
the contribution from the secondary Reggeons at LHC being presumably negligible. A 
comprehensive introduction to high-energy diffraction can be found in Ref. [8]. 

LHC will be the first accelerator where the relative contribution from secondary (sub- 
leading) trajectories (R) will be negligible, i.e. smaller than the experimental errors. The 
ratio R/IP, apart from kinematics, depends essentially on the difference of the relevant 
intercepts, however the above statement holds even for the most conservative (i.e. large 
R/IP) ratio, see e.g. Ref. [9]), decreasing with |t| since the Pomeron slope is smaller than 
that of the sub-leading Reggeons. 

Parametrization of the Pomeron is far from being unique. According to Ref. [10], there 
is only one Pomeron in the nature (although its form is not necessarily simple, see [11-16]). 
The data on deep inelastic scattering from HERA provoked discussions on the existence 
of an alternative, "hard" or "QCD Pomeron" [17-20], needed for the confirmation of both 
perturbative QCD calculations and of the "hard" and "semi-hard" diffractive physics. It 
should be remembered that the properties (parameters, etc.) of the Pomeron in hadronic 
collisions (ISR, SPS, Tevatron and LHC) can be determined with a precision and reliability 
much higher than that in ep collisions. 

The interface and/or transition between soft and hard dynamics is a key issue of the 
strong interaction theory. In elastic scattering at the LHC it is expected to occur in a 
smooth way, within the reach of the forthcoming LHC experiments. Roughly speaking, this 
region will be characterized by a transition from an exponential in t, through an exponen- 
tial in J\t\ ("Orear regime"), fall-off of the differential cross section to a power behavior. 



^It would be instructive and amusing to make a comparative compilation of earlier predictions already 
checked and thus confirmed or rejected by the existing data! 



manifesting hard scattering between point-like constituents of the nucleons. The transition 
region is not so simple because of the different unitarization and rescattering procedures 
used in the models. Moreover, non-linear trajectories (that mimic hard scattering), non- 
perturbative contributions and higher order perturbative QCD corrections complicate the 
issue. In particular, quark model and QCD calculations, apart from a power behavior in t, 
indicate [21] the onset of an s— independent regime in the differential cross section typical 
of the transition from soft to hard collisions. The complexity of this transition is connected 
with the deconfinement of quarks and gluons in nucleons. We argue in this paper that this 
transition is expected in the region (5 < — t < 15) GeV^, well within reach of the LHC 
measurements. 

We propose to investigate LHC effects that can be predicted qualitatively and which 
can be measured to give definite quantitative answers concerning the nature of the strong 
interaction at large distances. The LHC will rule out many of the existing model predictions 
and thus narrow the class of viable theoretical approaches. 

In Sees. 2 and 3 the essential features of the experimental program concerning forward 
physics at the LHC are described. 

In Sec. 4 we discuss two types of interesting irregularities observed in t within the 
exponential cone. One is the so-called "break" ^ of the cone (change of its slope) near — t = 0.1 
GeV^ . The second one concerns the possibihty of tiny oscillations superimposed over the 
cone. 

In Sec. 5 the "dip-bump" structure is analyzed. This is a sensitive region that will 
help to understand the nature of the high-energy diffraction and, eventually, to reveal the 
Odderon, whose role (see Ref. [22] and references therein) is often exaggerated, dramatized 
and confused. It is maintained (see e.g. Refs. [10, 14]) that the Odderon should exist 
simply because nothing forbids its existence. It is not known how large the Odderon is or 
how its contribution varies with s and t. Specific theorems and predictions concerning the 
Odderon can be found in Refs. [23-33]. Until now it has been observed directly in a single 
experiment [34, 35] only, and therefore needs to be confirmed experimentally. 

The region beyond the dip-bump structure, considered in Sec. 6, along the expected 
second cone, —t>l GeV^, may be indicative of the transition from soft to hard physics. 
It manifests transition from an exponential to a power decrease in t, with a possible slow- 
down of the energy dependence [21]. Details of this effect will measure the nonlinearity of 
the Pomeron trajectory at large enough \t\. One should note that the analytic 5— matrix 
theory, perturbative QCD and the data require that Regge trajectories be nonlinear, complex 
functions [36,37] (for more details see Refs. [10,38]). 

At the LHC energies a new phenomenon, namely the onset of the Black Disk Limit (BDL) 
may come into play, changing the t-dependence of the slope B{s,t). This phenomenon is 
discussed in details in Sec. 7. 

The Pomeron trajectory has threshold singularities, the lowest one being due to the two- 
pion exchange, required by the t— channel unitarity [36,37,39-42], as shown in Fig. 1. There 



^Actually, the "break" is an approximation to a smooth curvature. 
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Figure 1: Two-pion loop, required by t— channel unitarity, in the pomeron exchange, 
is a constrain [36,43-46], following from the t— channel unitarity, by which 

53a(^)~(^-^of°(W/^ t^to, (i) 

where to is the lightest threshold in the given channel. For the Pomeron trajectory the 
lightest threshold is to = 4m^, as shown in Fig. 1, and the trajectory near the threshold can 
be approximated by a square root: 

ait) ~ ^Aml - t. (2) 

The observed nearly linear behaviour of the trajectory is promoted by higher, additive thresh- 
olds (see Ref. [47,48] and references therein). 

This threshold singularity appears in different forms in various models, see Sec. 4. It is 
also noted that, irrespective of the specific form of the trajectory or scattering amplitude, the 
presence of the above-mentioned threshold singularity in t results in an exponential asymp- 
totic decrease of the impact parameter amplitude, with important physical consequences 
(the so-called nucleon atmosphere, or its clouding). 

The important role of nonlinear trajectories and their observable consequences were first 
studied in Refs. [37,49-52]. Independently, they were also developed by the Kiev group (see 
Ref. [53-55] and references therein), and more recently in a series of papers (Ref. [56] and 
references therein). 

Asymptotically, the trajectories are logarithmic. This asymptotic behaviour follows from 
the compatibility of the Regge behavior with the quark counting rules [57] as well as from 
the solution of the BFKL equation [17-20]. A simple parametrization combining the linear 
behaviour at small \t\ with its logarithmic asymptotic is [10,38,58] 

a(t) = ao -7ln(l -/5it). (3) 

Such a trajectory, being nearly linear at small reproduces the forward cone of the dif- 
ferential cross section, while its logarithmic asymptotic provides for the wide-angle scaling 



behavior [57-59]. Eqs. (2) and (3) can be combined in the form [10,14] 

a{t) = ao-i ln(l + /^aVV^) , (4) 

where jSi, /?2 and 7 in Eqs. (3) and (4) are parameters whose numerical values can be 
associated with their physical meaning (see [60]). 

In a limited range, especially at small and intermediate values of their argument, linear 
trajectories may be a reasonable approximation to their otherwise complex form. 



2 Measurement Strategy at the LHC 

The elastic proton-proton interactions are measured, and triggered, by the leading proton 
detectors, Roman Pots, placed symmetrically on both sides of the CMS experiment at ±147 
and ±220 meters from the Interaction Point (IP5) ^. To measure protons at small scattering 
angles, the detectors must be moved close to the primary LHC beam in vertical direction ^. 
The "nominal" TOTEM beam optics set-up has high /3* ( /3* ~ 1540 m) and no crossing 
angle at the IP for optimizing the acceptance and accuracy at small values of the four- 
momentum transfer squared down to — tmm ~ 2 ■ 10~^ GeV^. 

The —t distribution of the scattered protons, dNei/dt, is extrapolated to —t = 0, where it 
is related to the total proton-proton cross section by the Optical theorem. With the special 
optics of P* = 90 m, compatible with the LHC injection optics, a first quick measurement 
of the elastic cross section daei/dt could be made; extrapolation to the Optical point will be 
made with an accuracy of a few percent. 

Runs with different LHC optics set-ups, such as (3* = 90 m, the "injection" optics 
= 11 m), stage-1 "pilot" run optics { (3* = 2 m) and the standard LHC optics 
{(3* = 0.55 m) will allow measurements up to — t ~ 10 15 GeV^. Runs with a reduced 
center-of-mass energy will allow an analysis of energy dependence, comparisons with the 
Tevatron results and a precise measurement of the p parameter. 



2.1 Elastic proton signature at the LHC 

Transverse position of an elastically scattered proton with the momentum loss C, = Ap/p at a 
distance s from the interaction point (IP): ( x{s), y{s)), is given by the initial coordinates at 
the IP, {x*{s = 0), y*{s = 0)), scattering angle, 9* y, the effective length L^^J , magnification, 
Vx^y, and dispersion, D, as (see Ref. [61] and references therein) 

•^The detector locations at ±145 and ±149 meters from IPS are here referred as the "±147" meter location 
and the ones at ±218 and ±222 meters from IP5 as the "±220" meter location. Initially the ±220 m location 
will be instrumented. 

^The closest approach to the beam is of the order of a few mm's (lOcr + 0.5 mm) and depends on the 
measurement location and used beam optics scheme (see Ref. [61] ). 



By using the transfer matrix, T(s), the initial coordinates of an elastically scattered proton at 
the IP can be mapped into a detector location at point s along the machine. By measuring the 
transverse position and scattering angle at point s, the initial coordinates can be determined. 
Since angles 6* y are very small, they have to be measured by combining the symmetrically 
located detector stations on both sides of the IP: a combined measurement, using the left and 
right arms of the leading proton detectors, yields an accurate measurement of the coUinear 
pairs of elastically scattered protons (Fig. 2). 

From Eq. (5) 



xr - XL 



Vr - Vl 



(6) 



The elastically scattered protons have to be measured close to the primary LHC beam; 
down to angles of 5 to 10 /xrad with respect to the beam direction. For this, special beam 
optics conditions with reduced beam divergence at the interaction point (IP) - and sufficiently 
large displacement of the scattered protons at the detector locations (Fig. 2) - are required. 
Thorough studies of optimal LHC beam conditions have led to the nominal "TOTEM" beam 
optics with (3* ~ 1540 m ^ for the measurement of elastic scattering and soft diffraction 
Ref. [61]. 
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Figure 2: The TOTEM lay-out of leading proton detectors (Roman Pots). The detector 
locations at ±147 m (RPl) and at ±220 m (RP3) are shown [61]. The location of RP2 is 
not part of the present design. 

In Fig. 3, the acceptance of elastic protons is shown as a function of —t for three different 
run scenarios. An acceptance down to ~ 5 ■ 10~^ GeV^ is achieved with the nominal 
LHC beam emittance; with improved beam emittance an acceptance down to — tmin = 
2 ■ 10"^ GeV^ is obtained [62]. Contrary to the nominal TOTEM optics, the /3* = 90 m 
optics uses the standard LHC injection optics and could be realised during the run-in phase 
of the machine. The acceptance in —t reaches —tmin ~ 3.0 ■ 10~^ GeV^. 

Due to the geometric constraints imposed by the LHC vacuum pipe and beam screen, 
elastic protons with —t values in excess of 0.5 GeV^ cannot be measured with the nominal 
high-/?* optics conditions. Also, since at high —t values daei/dt oc 1/t®, higher luminosities 
would be desired. The "standard" LHC optics scenario includes the "injection" optics with 



^There is significant uncertainty in determining the exact value of /?*, especially at large values of /?*, 
and the number given should be taken as an approximate figure, only. 
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Figure 3: —t acceptance of elastic protons detected at ±220 m from the IPS for the 'nominal' 
(/?* = 1540 m) and for two 'custom' run options [jS* ~ 90 m.) [62]. 

/9* = 11 m and a "pilot" stage optics with (3* = 2 m that allows elastic protons with higher 
—t values to be accessed (Figure 3). With the 2 m optics, a reasonable event statistics up 
to -t = 10 GeV^ (15 GeV^) and resolution a{\t\)/\t\ ^ 2% is achieved. 

With the high-/?* TOTEM optics, so-called parallel-to-point focussing conditions are 
achieved (at ±220 m in both horizontal and vertical planes) and the elastic proton mea- 
surement becomes independent of the location of the primary interaction vertex. This will 
reduce systematic uncertainty in measuring the four-momentum-transfer squared of the scat- 
tered proton. 

The angular beam divergence dominates the uncertainty in measuring -t, and together 
with the detector resolution, accounts for basically all of the uncertainty in —t. For small 
values of —t {—t = 0.01 GeV^) the error in —t, At/t, is ~ 3% in case the forward-backward 
pair of leading proton measurement stations is used and ~ 5% if only one of the two leading 
proton spectrometry "arms" is used. 

With lower center-of-mass energies, the acceptance in —t improves and, at the energies of 
^/s = 2 TeV, —t values of one order of magnitude lower than at the nominal LHC energy are 



reached. This would allow comparisons with the Tevatron results and the Coulomb-hadronic 
interference region to be probed. 



2.2 Extrapolation to the optical point 

For the total pp cross section measurement, extrapolation of the elastic scattering t distribu- 
tion, dNei/dt, to t = is required. The relative statistical uncertainty of the extrapolation 
is estimated to be 0.1% based on short periods of data taking with the nominal TOTEM 
optics and luminosity of 10^^ cm~^s~^ [61]. 

A thorough study of the systematical effects in the extrapolation process was carried out 
[61,62] and concluded that a precision better than 0.5% should be achieved in extrapolating 
the elastic cross section to the optical point^. 

In case Coulomb scattering is not accounted for in the extrapolation process, a shift of 
{daei/dt)\t=Q of (1 -7- 2)% could occur. Due to Coulomb scattering contribution, the slope B 
in da^i/dt ~ exp{Bt) does not stay constant at low —t as assumed in a usual extrapolation. 
The uncertainty due to this effect might represent an uncertainty of several parts in 10"'^ to 
the extrapolation. 



2.3 Measuring p 

For measuring the p-parameter. 



one fits the elastic differential cross section: 

L dt dt t2 iQQ^ ' ^ ) 

where the three terms are due to Coulomb scattering, Coulomb-hadronic interference, and 
hadronic interactions. L is the integrated luminosity, a the fine structure constant, A0 the 
relative Coulomb-hadronic phase, given as 

A0 = -ln(5(s)|t|/2) -0.577 (9) 

and Fi{t) is the nucleon electromagnetic form factor, which is usually parameterized as 

In the least squares fit procedure, the following two equations are also used: 

2 1 167r dNei, 



L(l+p2) dt 



|t=0) 



:iii 



^For a discussion on model dependences in extrapolation to the Optical point see Ref. 



(12) 



Eq. (11) is a direct consequence of the Optical theorem. N^i is the total number of events 
obtained by integrating the dNei/dt distribution within the —t region where hadronic inter- 
actions dominate, and extrapolated to — t = and to t ^ oo by using the form exp(— i?|t|). 
Ninei is the total number of inelastic events. Note that Eqs. (11) and (12) allow the lumi- 
nosity L to be expressed in terms of atot and p. Then dNg,i/dt in Eq. (8) can be expressed 
in terms of just three unknowns: atoti B and p. In the fit procedure, the same data on 
dNei/dt, together with the total number of inelastic events Ninei recorded during the same 
experimental data taking runs, are used as inputs. A least-squared analysis for at^t, B and 
p in Eq. (8) is done by using all the collected input data. 

The evaluation of systematic errors due to the uncertainty in beam emittance, vertex 
positions and spread, beam transport and incoming beam angles is based on Monte Carlo 
and machine simulations. These simulations use the geometry of the experimental set-up 
and efficiency of the detectors as input. 

2.4 Elastic scattering run scenarios 

For elastic scattering, three run scenarios are considered (Table 1): 

1. Nominal TOTEM optics for (low —t) elastic scattering, (3* ~ 1500 m, 

2. An early medium- (3* optics, with (3* = 90 m, and 

3. Optics for large —t elastic scattering, f3* = 2 — 0.55 m. 

The event rate per bunch crossing is calculated as (for symbols see Table 1) 



where aei = elastic cross section, L = luminosity, / frequency. A'';, = no. of bunches; factor 
{Ns/Nh ~ 1) accounts for the empty buckets. 



2.5 Total pp cross section measurement strategy 

The total proton-proton cross section is measured - in a luminosity independent way - by 
using the Optical theorem. By extrapolating the elastic rate down to the Optical point, t = 0, 
and by recording the elastic and inelastic event rates, the total cross section is measured with 
an over-all accuracy better than 1%. 




(13) 



Scenario 


1 


2 


3 




low \t\ elastic, 


low \t\ elastic. 


large \t\ elastic. 


Physics: 


1%), 

MB, soft diffr. 


(Ttot(@~ 5%), 

MB, soft diffr. 


hard diffraction 


P*[m] 


~ 1500 


90 


2 ^ 0.5 


N of bunches 


43 ^ 156 


156 


936 ^ 2808 


Jjunch spacing [nsj 


2025 ^ 525 


525 


25 


N of part, per bunch 


(0.6^ 1.15) ■ 10" 


1.15 ■ 10" 


1.15 ■ 10" 


Half crossing angle [firad] 








92 


Transv.norm.emitt. e„[/imra(i] 


1 


3.75 


3.75 


RMS beam size at IP [fim] 


450 


213 


32 


RMS beam diverg, at IP [firad] 
Peak Luminosity [cm~^s~^] 


0.3 


2.3 


16 


1028 ^ 2 . 10^9 


3 • lO^o 


1033 



Table 1: Three different LHC run scenarios foreseen for elastic scattering measurements at 
the LHC. 



For the total cross section measurement, the "nominal" TOTEM beam optics (/?* ~ 1500 
m) with several short runs is used. During the initial LHC running, the run scenarios with 
/3* ~ 90 m is planned to be used for a total cross section measurement with an accuracy of 
about 5%. 

For measuring the total cross section, inelastic scattering needs to be studied within large 
Et- The aim of the forward physics initiatives at the LHC is to complement the base line 
ATLAS (ALFA), CMS (TOTEM) and ALICE designs with forward detector systems (see 
Refs. [64-68] ). Besides restricted detector acceptance, inadequate theoretical understanding 
of the forward physics phenomena poses a serious systematic uncertainty for the base line 
experiments in need of precise luminosity measurement. As an example, the single diffrac- 
tive cross section for the low diffractive mass region, M* < 5 GeV, could amount to 25% 
of the over-all agd and cause a major systematic uncertainty in the total pp cross section 
measurement. 

The base line LHC experiments define a "minimum bias" event category that must be 
suppressed due to the limitations in recording pp interactions in excess of the rate (P* level 
trigger band width) which can be recorded, as the cost of triggering on "interesting" large 
Et events. The main task of equipping the forward region of a main stream LHC experi- 
ment is to complement the physics reach by including the events that are not selected by the 
"minimum bias" event trigger. The TOTEM experiment together with CMS forward detec- 
tors (CASTOR, ZDC, and the proposed FSC and FP420 systems), represent the necessary 
complement for selecting an unbiased (!) sample of "minimum bias" events required for the 
analysis goals stated in the CMS Physics TDR. 

In general, the soft particles in the non-diffractive event category (nd) will end up at 
central rapidities, while the relatively few energetic ones are expected to end up at small 
angles, to be recorded by the forward calorimetry and spectrometers. 

The longitudinal, z, coordinate of the vertex is defined by first determining the distance 
of closest approach to the nominal z axis for each track, Ztrack, (at least two tracks are 
required) and by then calculating the mean: z^x =< Ztrack > ■ By using this simple method. 



the z coordinate of the primary interaction vertex can be determined with a resolution of 
Ovx ~ 5 cm. The beam related events are identified by requesting that their reconstructed 
z^x value is within ±50 cm of the nominal IP. This selection is found to be more than 96% 
efficient in choosing the beam related events. The remaining events have diffractively excited 
systems with relatively small masses, M* < 10 GeV, in which all the charged tracks escape 
detection in TOTEM Tl and T2 spectrometers (Figure 4). 
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Figure 4: Left panel: Ratio of detected single diffractive (sd) events as a function of diffractive 
mass. Right panel: uncorrected Monte Carlo simulated (unshaded) and acceptance corrected 
(shaded) sd events with diffractive masses above 10 GeV [61]. 

For obtaining the over-all inelastic rate, the missing coverage''' above i] > 7 was estimated 
by extrapolation. Fig. 4 (right panel) shows the simulated 1/M*^ distributions for the single 
diffractive (sd) events before and after acceptance correction. A linear fit is used to correct 
for the unseen part of spectrum. In the case of sd events, a correction of 4%, corresponding 
to about 0.6 mb, was quoted above the detected fraction of this category of events. A similar 
analysis gives a correction factor of 0.1 mb for double diffractive and 0.02 mb for the central 
diffractive events [61]. 

Unfortunately diffractive cross sections are poorly known and theoretical understanding 
of both small and large mass single diffraction is lacking. Low mass single diffraction (M* < 
5 GeV) could represent 25% of the total diffractive cross section (for double diffractive 
{dd) events the uncertainties are even more severe). Moreover, in the light of recent CDF 
measurements, soft central diffractive cross sections are likely to be seriously overestimated 
in current Monte Carlo models. 

With the vertex constraint, a substantial part of the beam-gas events is rejected. The 
simulation studies [61] show that by using the vertex constraint, the beam-gas interactions 
closer than ±5 m represent only ~ 3% of the selected sample of events. Since the beam-gas 
event rate is estimated to be 2 Hz/m at 1/20 of the nominal beam current, the trigger rate 
due to this source of background can be safely neglected. In addition, the Roman Pot based 
leading proton trigger (e ~ 90%) can be used to further reduce the background rate. 



^The rapidity coverage could be extended further by simple Forward Shower Counters (FSCs) placed at 
60 to 140 meters from IP5. 



2.6 atot measurement 



The LHC (TOTEM) measurements of the total pp cross section is luminosity independent 
and based on the Optical theorem: 

16 TT {dN/dt)t=0 

(l + p2)iVe, + iV_,' ^'^^ 

where A'e; and Ninei are the elastic and inelastic event rates, {dNei/dt)\t=o is the elastic cross 
section extrapolated to the Optical point. 

The relative error in atot, neglecting uncertainty in p , is then: 



datot\^ 



( d{dN,i/dt)t=o ^ ^ / d{N,i + Nir^el] 
(ytot J [ idN,i/dt)t=o J [ (iVd + N.nel) 



(15) 



The uncertainty in inelastic cross section is estimated to be less than 1 mb (see Ref. [61]). 
However, little experimental data on small mass diffraction exists, and the uncertainty in 
inelastic diffraction alone could amount to several millibarns. Together with the uncertainty 
of about 0.5% of the extrapolated value of {daei/dt)\t=o, results in an over-all error of 



(^tot 



» 1%. (16) 



The uncertainty in the value of the p-parameter could also have an important contribution 
to the atot measurement (see Chapter 3 and Ref. [69]), and it could be reduced by a direct 
measurement. 



2.7 Luminosity measurement and monitoring 

Luminosity measurement 

The luminosity relates the cross section, ai, of a given process i to the corresponding 
event rate Ni by 

i=^. (17) 

It is not trivial to find a process with a well defined - and precisely calculable - cross section 
combined with a prominent event signature. The cross section should be large enough for 
monitoring the luminosity as a function of time, e.g. during a fill or when investigating bunch- 
to-bunch variations. By using several complementary aprroaches, systematic uncertainties 
can be brought under control [63]. 

A simultaneous measurement of the elastic and inelastic event rates can be used to define 
the luminosity as 

16 vr {dNjdt)\t=o)' ^ ' 



On-line luminosity monitoring 

LHC is the first hadron collider where, due to the large c.m.s energy and high luminos- 
ity, a significant number of inelastic interactions (on average 35 interactions at the design 
luminosity) are expected to take place per bunch crossing. The traditional technique of 
monitoring luminosity by requesting a coincidence of two counters at small angles on both 
sides of the IP is not sufficient at the LHC. At the design luminosity, the coincidence rate 
observed will not be proportional to the number of events, but rather to the number of hunch 
crossings. In the case of no segmentation of the luminosity monitors, the probability to get 
a coincidence will be close to 100%. 

Highly segmented forward detectors in T2 and/or FSC could be used as luminosity mon- 
itors. The rate used to monitor the luminosity, could be defined by using the double-arm 
coincidence rate between a pair of left-right detector segments. The segmentation reduces 
the counting rate significantly below the bunch crossing frequency and, therefore, becomes 
proportional to the luminosity. A coincidence signal between a pair of left-right detector 
segments is unlikely in case of separate overlapping pp collisions. The technique would help 
to suppress beam related backgrounds. In principle, beam related background may be rec- 
ognized by a time stamp given by a forward detector. During the beam cross-over at the IP, 
no other bunches should pass through the luminosity monitor location. Secondary particles 
from beam-gas and beam-wall interactions are traveling in-time with the bunch. In practice, 
the time stamping is challenging due to the high bunch crossing frequency (40 MHz). By 
requiring simultaneous left-right signals contributions from the beam related backgrounds 
are practically eliminated. 

The calibration of the luminosity monitors can be performed during the dedicated high- 
13* runs at lower luminosities, where the luminosity is precisely determined together with 
tot. Alternatively, once the total inelastic cross section is precisely measured during the 
high- [5* runs, together with the elastic and total cross sections, the inelastic events could be 
used to re-calibrate the monitor during the \ow-(3* running. This becomes important when 
any significant changes to the detector lay-out are made, e.g. when the outer detectors are 
dismantled after the first year of running. 

3 Forward physics at the LHC 

Forward, or soft physics, roughly speaking, is the synonym of diffraction - elastic and inelas- 
tic, the role of the latter increasing with increasing scattering angle (or momentum transfer). 

Forward physics will play an important role in early runs of the LHC for at least two 
reasons. One is that measurements of basic quantities, such as the total cross section (Jtot{s), 
the ratio p{s) of the real part to the imaginary one of the forward scattering amplitude, 
the local slope B{s,t), of the differential cross section, etc. are of fundamental importance 
for the calibration/normalization of the beam and detectors and this task goes beyond the 
problems of understanding the nature of diffraction. Secondly, apart from classical studies of 
diffraction (Pomeron), the diffractive medium (gluons) may also favour central production 
of the Higgs boson. With increasing luminosity, the experiments at LHC will gradually shift 



towards measurements of rare events in the non-forward direction. 

At the LHC, three collaborations, namely TOTEM/CMS, ALICE and ATLAS, are 
preparing to measure elastic, inelastic and the total pp cross section * at the expected en- 
ergy a/s = 14 TeV. The complementarity of the expected performances will ensure optimal 
reliability of the results. The total cross section is claimed to be measured within 1% of 
precision. The precision of the extrapolation to the optical point has been analyzed in [70] ^ 

One should note that the predictions based on model extrapolations for afo^ at -^i = 14 
TeV have a wide range. For example, Landshoff predicts [1,2] afot (14 TeV) = (125±25) mb. 

TOTEM, CMS and FP420 collaborations are combining their efforts to cover a phase 
space (see Fig. 2), where the geometric acceptance of detectors is shown) exceeding that in 
any of the preceding collider experiment [67, 71, 72]. 

The TOTEM experiment, in particular, aims at measuring 1) the total proton-proton 
cross-section with a relative precision of about 1%; 2) the elastic proton-proton scattering 
between 10~^ GeV^ < —t ^ (p0)^ < 10 GeV^, p and 6 being respectively the proton 
momentun and scattering angle in the c.m.s. Furthermore, by combining TOTEM with 
CMS and its forward additions, the CASTOR and ZDC calorimeters, a full forward physics 
program is planned. 

The TOTEM experiment will measure pp elastic scattering to about — t = 0.15 GeV^ and 
will provide crucial new understanding of the phenomena in the high —t elastic scattering 
regime. To discriminate between different models, it is important to measure elastic scat- 
tering in the widest possible kinematic region. In most of the previous measurements of the 
total cross section (at ISR, SPS, Tevatron, RHIC) the value of the parameter p was imported 
either from another experiment or from model calculations. As argued in Refs. [72,73], known 
variations of the value of p have a negligible effect on the resulting value of atot- This is in 
disagreement with Ref. [74], where it was argued that a small variation of p may affect the 
resulting atot significantly. If this is true, a simultaneous fit to all inputs {o-toti P, and C) 
from a single experiment is needed. Publication of direct and unbiased data on sd — dN / dt 
is highly welcome. 

There are two possible kinds of luminosity measurements: one yields an absolute value 
which serves as a point of reference, the other one gives a relative value as a function of time. 
The latter measurement will be performed by ATLAS using a special detector called LUCID 
(Luminosity measurements Using Cherenkov Integrating Detectors). The idea of the LUCID 
detector is explained e.g. in Ref. [75]. It makes possible measurements in the very forward 
region, directly related to the instantaneous luminosity. This detector will enable ATLAS 
to obtain a linear relationship between luminosity and the number of tracks counted in the 
detector, directly related to the luminosity (see Table 2). 

Measurements in the Coulomb-nuclear interference (CNI) region, 

~t ~ 6.5 ■ 10"^ GeV^ Q^in ~ 3.5 prad, 

(note that at SPS Qmin ~ 120 prad) will be be used to extract the value of the parameter p 

^ALICE will use TOTEM's measurements of the total cross seetion. 

^There are some reservations about the achievable precision due to the uncertain cross section of low-mass 
diffractive scattering that could amount to several millibarns. 
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101.14 mb 


0.9% 


B 


18 GeV-2 


17.93 GeV-2 


0.3% 


P 


0.15 


0.143 


4.3 % 



Table 2: Result of a fit (see Ref. [76]) to simulated dN/dt data corresponding to ~ 1 week 
(10 M events) of running. C = 10^^cm~^s~^. 
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Table 3: Range of the measured values of atot above the ISR energies 
from Eq. (8). 

The extracted value of p may be affected by at least two phenomena. One is connected 
with the well known corrections [77] to the Coulomb-hadron phase [78] and the other one 
with the non-exponential behavior of the diffraction cone, known as the "fine structure of the 
Pomeron trajectory" [23-25,39-42]. Below, in Sec. 4, we shall come back to this point. The 
program of studying this "fine structure of the Pomeron trajectory" is among the priorities 
of the ATLAS collaboration [73,75]. 

The ALICE experiment is designed as a general purpose experiment with a central barrel 
covering the pseudorapidity range —0.9 < r] < 0.9 and a muon spectrometer covering the 
range —4.0 < t] < —2.5 at luminosities C = 5 ■ 10^°cm~^s~^ and C = 10^"^ cm~^s~^ in 
pp and PbPb collisions, respectively, as well as an asymmetric system pPb at a luminosity 
C = 10^^ cm-h-^ (see, also, Ref. [79-81]). 

Moreover, the experimental program of ALICE to large extent will be oriented to inelastic 
reactions, e.g. by studying their dependence on the width of the rapidity gap. 

It should be noted that all these values are approximations since some of them were taken 
at t = as p(0), dNei/dt \t=o, or integrated over all the region of angles as Nei, Ninei, i.e. all 
these values were obtained under some theoretical assumptions. 

Tables 3 and 4 [82] show the range of the measured and the expected values of atot above 
the ISR energies, with the divergence of the existing data on atot at the same energy quoted 
in Table 2. One can see that the extraction of atot from the differential cross section dN/dt is 
a complicated problem and that the use of different models can chance the predicted value of 
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Table 4: Predictions of the parameters of elastic scattering amplitude at {^/s = 14 TeV, 
t = 0) 



the total cross sections at the LHC energies from 80 mb to 230 mb (see Table 3). Therefore, 
one cannot use the theoretical predictions of, say atot or p, to extract other observables from 
the experimental data on dN/dt. With the exception of the UA4 and UA4/2 Collaborations, 
the numerical data on dN/ dt of other collaborations were not published, excluding any cross- 
check or improvement of these results. We hope that the future LHC data on dN/ dt will be 
published. 

In Ref. [75] the significant correlations between the values of p and atot are also visible. 
The results of the fits of the simulated LHC experimental data in the framework of the 
non-exponential model of the hadron scattering amplitude presented at "EDS-07" [82] at 
a/s equal to 2 TeV and 14 TeV show large errors in the determination of atot and p(0) (see 
Tables 3 and 4). 

The (weak) dependence of atot on p, if p is very small, and it comes from the coefficient 
l/(l-|-p^) in front of da/dt (see Eq. (14)), while the strong dependence of the normalization of 
dN/dt, the values of atot and p(0) comes from the extraction of the Coulombic and Coulomb- 
hadron interference terms from dN/dt data to obtain {daf.i/dt)t=o- The Coulomb-hadron 
interference term is proportional to [p{s,t) — a^cn{s,t)]. At very small t, the Coulomb- 
hadron phase is also important. Corrections to the Bethe formula [94] were calculated in 
Refs. [95-98]. A detailed analysis of the role of this term in the behaviour of the differential 
cross sections was carried out in Ref. [92]. To make the analysis complete, the Odderon 
contribution [23-25] as well as the nearby threshold singularity at t = 4m^ [39], should be 
also taken into account. 

To minimize the errors, in some experiments the value of atot or p were fixed from other 
measurements. This was done, for example, by the UA4/2 Callaboration, which extracted 
p(0) by using atot from the result of the UA4 Collaboration [atot = 61.9 mb). However, the 
value of atot obtained by the Collaboration UA4/2 turned out to be atot = 63.0. With such a 
value of atot the resulting p(0) becomes larger. This situation was analysed in Refs. [69,99]. 
with the results shown in Table 5. 





p = 540 GeV, 


0.000875 < \t\ < 0.12 GeV2) 


n 


experiment exp. analysis 


mod. anal. I [99] mod. anal. II [69] 


1 


UA4 0.24 ± 0.02 


0.19 ±0.03 


2 


UA4/2 0.135 ±0.015 


0.17±0.02 



Table 5: Comparison of experimental values of the ratio p of the real to imaginary part 
of the scattering amplitude, obtained in the UA4 and UA4/2 experiments, with theoretical 
values. 

Diffraction dissociation, in particular the low mass one, is among priorities of the first 
experiments at the LHC It should be remembered that diffraction dissociation - both single 
and double - and elastic scattering have much in common. Similarities are expected in 
the shape of the diffraction cone, with its fine structure (the "break" and oscillations), 
open for observation, see Ref. [100], and the dip-bump structure (until now not seen in 
diffraction dissociation). Hence studies of elastic scattering and diffraction dissociation are 
complementary. A recent overview of the ALICE detector and trigger strategy for diffractive 
and electromagnetic processes at the LHC can be found in Ref. [79-81]. 



4 The forward cone 



4.1 The "break" at very small \t\ 

An essential part of the future TOTEM±CMS and ATLAS experiments is connected with 
the measurement of the elastic differential cross section at small momentum transfer, with 
the purpose of extracting from these data the values of the total cross section. An important 
point of this procedure is the simultaneous measurement of four quantities: the luminosity 
L (or the normalization coefficient), the total cross section atot, the slope B{s,t), defined as 

and the ratio p. The last two quantities depend on their, a priori unknown, t-dependence. 
Consequences of this complexity are the contradictions of the obtained values of p{s) in differ- 
ent experiments in the energy range of a/s = 5-^20 GeV [99], in particular the large difference 
in the values of p as measured in the UA4 and UA4/2 experiments at ^/s = 540 GeV as 
well as the discrepancy between the values of atot at -\A = 1-8 TeV. This discrepancy results 
from the correlations between simultaneous measurements of several unknown values in a 
single experiment instead of fixing some values from the phenomenological analysis or from 
other experiments. 

This situation was demonstrated for example in Ref. [99]. Since we do not know the 
true t-dependence of p and of B, one cannot use the same constant values for different 
small intervals of t. Therefore, simultaneous measurements of all unknown observables in 
a single experiment are needed, rather than their substitution by values taken either from 
phenomenological models or alternative experiments, as illustrated e.g. in Refs. [69,99], 
where compatible values of p, different from those experimental were obtained. 



Given the uncertainty in the s and t dependence of the observables, the standard pro- 
cedure is to suppose their smooth and monotonic behavior from low to super-high energies. 
Along these lines, it is usually proposed to measure atot taking the value of p from the results 
of theoretical analysis (for example, that of the COMPETE Collaboration [86]), although, as 
argued in Refs. [3-6] , any consistent analysis should include a simultaneous fit to expression 



In the ISR energy region, the pp diffraction cone changes its slope B[s,t) near t = 
—0.1 GeV^ by about 2 units of GeV~^. In Ref. [39] this phenomenon was interpreted as 
the manifestation of t— channel unitarity (a two-pion loop, see Fig. 1) and, in term of the 
Pomeron exchange, it was modeled by the inclusion of a relevant threshold singularity in the 
trajectory^^: 



where ai ~ 0.1. It will be interesting to see whether this effect will persist at the LHC 
and thus confirm the idea of Ref. [39] . 

As shown in Ref. [100], the "break", or the fine structure of the Pomeron, in principle, 
can be seen directly, unbiased by electromagnetic interactions, in proton-neutron scattering. 

4.2 Oscillations 

An important phenomenon that should be scrutinized at the LHC is the possible appear- 
ance of small-period oscillations over the smooth exponential cone. Possible deviations from 
a simple exponential behaviour of the hadron-hadron scattering amplitude were discussed 
in the literature long ago. For example, it was shown in Ref. [101] that peripheral contribu- 
tions from inelastic diffraction result in large- and small-period oscillations in the momentum 
transfer. Among the attempts to verify these oscillations the first one was done in an experi- 
ment at Serpukhov, where oscillations in pp scattering were detected [102] and then discussed 
in Ref. [103]. An alternative view, interpreting these oscillations as an artifact, connected 
with the t— dependence of the slope, was put forward in Ref. [104]. In Ref. [105] the statis- 
tical nature of the possible oscillations in the ISR data were analyzed in the framework of 
the Dubna Dynamical (DD) model (for the DD model see Subsection 5.4). The results are 
shown in Fig. 5. 

The effect was first observed at the ISR and subsequently was discussed in Ref. [106,107]. 
Although it has not yet been confirmed unambiguously, it continues to attract attention (see 
Ref. [106,107] and references therein). The analysis of the slope [106,107] shows possible 
oscillations in the UA4/2 data. The super-fine structure (oscillations) of the cone may be 
related to residual, long-range interaction between nucleons [108] or the action of a potential 
of rigid hadronic strings [109]. 

A new method of the analysis of the experimental data of UA4/2 experiment was pro- 
posed in Ref. [110]. The method is based on the comparison of two statistically independent 



of Eq. (8). 




(20) 



This threshold smgularity can be included also in the logarithmic trajectory as in Eq. (4). 




Figure 5: Ratio of the experimental differential cross section to theoretic one calculated from 
the DD model [105] in the ISR energy region. 



sets, Ref. [111]. If we have two statistically independent sets and x^^ of values of the 
quantity X distributed around a certain value, say A, with the standard error equal to 1, we 
can try to find the difference between x^^ and x"^^ . For that we can compare the arithmetic 
mean of these choices: 

AX = (x'l + X2 + ■..x'^i)/ni - {xl +xl + ...xl2)/n2 = x'^^ - x^^. 

The standard deviation for this case will be 

fe= [iM + l/ns]^/^ 

If we have the purely statistical noise the value AX/5x tends to zero. However if there 
is some additional signal, this value should differ from zero. When it is larger than 3, one 
can say that the difference between these two choices has a 99% probability of confirming 
the presence of the oscillations (for more details see Ref. [110]). 

The deviation Ai?j of each of the experimental data for the cross section from the corre- 
sponding theoretical values is measured in units of the experimental error 51^-^: 

ARi = [{da/dtiY'^P - {da/dtif'^ydT''. (21) 



By summing these ARi over all 99 experimental points of the UA4/2 experiment, the 
result should tend to zero as the statistical deviations are equally probable in both sides 



of the theoretical curve. However, if the theoretical curve does not precisely describe the 
experimental data, for example, the scattering amplitude deviates from the exponential 
behavior in the momenta transfer, the sum over Ai?j can differ slightly from zero, going 
beyond the value of a statistical error. To take into account this effect, we divide the whole 
interval of the momentum transfer into k equal pieces of size A such that kA > {qgg — qi), 
where = ^J\ti\, and sum ARi separately over the even and odd pieces. Thus, we get two 
sums and L"^^ for the rii even and n2 odd interval, respectively. For rii + n2 = k and 
|ni — 71.2! =0 1- 

ni 99 n2 99 

^up ^ Ai?i)|Ag(2j-l)<g,<Ag(2j); = '^C^ ^Ri)\Aq{2j)<q,<Aq{2j+l)- (22) 

j=l i j=l i=l 



In Ref. [110], where this method was applied to the data of the UA4/c Collaboration, 
Ag = 0.9085 ■ 10"^ GeV was used. 

Let us calculate the quantities L"^ and L"^"; the results are shown in Fig. 6 by the full 
and dash-dotted lines. It can be seen that in the range < \t\ < 0.1 GeV^ these quantities 
change drastically and in the range |t| > 0.1 GeV^ instead they vary slightly. It means that 
the amplitude of a possible periodic structure is decreasing with growing t. 




Figure 6: Left panel: Sums L"*' and L'^'^ calculated [69] with atot = 62.2 mb; B = 
15.5 GeV~^; p = 0.135 for go and Sq (full and dashed lines); and for go + Sq/2 ( dots and 
dots-dashed lines). Right panel: Same as on the left panel but with atot = 63.54 mb; B = 
15.485 Ge1/-2; p = 0.158. 



Note that this new method can be used to check the true determination of the parameters 
of the elastic scattering at small t. The two curves obtained have to be symmetric with 
respect to the line calculated by using the basic parameters {atot,P,B). 

A way to identify possible oscillations at the LHC is to use the method of overlapping 
bins, suggested by J. Kontros and Lengyel in Ref. [106]. The procedure consists in scanning 
the cone by overlapping bins in t, each containing a certain number of data points shifted 
by a small number of points n. Within each bin an exponential fit is applied [106, 107] (see 
Fig. 7). 




Figure 7: Left panel: t— oscillations calculated in Ref. [106] as the ratio of the experimental 
to the theoretical values at y/s = 19.4 GeV. Data are from Ref. [112]. Right panel: Same as 
on the left panel but for the data from Ref. [113]. 



The length of the bin should not be too large to oversimplify the parameterization and not 
too small in order to contain a reasonable number of points for each process. For example, 
a shift from bin to bin 6t = 0.01 GeV^ was found in Ref. [106] reasonable. 

Looking for oscillations in the cone at LHC by using the method of overlapping bins of 
Ref. [106, 107] is a promising program for future experiments. As already mentioned, the 
oscillations can be related to residual long-range forces between nucleons [108]. 



5 The Dip-Bump Region, t ^ — 1 GeV^ 

Before going into details, we would like to notice a model-independent regularity found in 
Ref. [114]. As shown in that paper, a correlation between the value of p and the depths of the 
minima of the diffraction differential cross section in proton-proton and proton-antiproton 
scattering exists. The ratio p for the latter changes sign in the energy region ^/s = 9.8 GeV. 
At this energy the differential cross section of pp-scattering in the dip-bump structure has 
its sharpest minimum, while p in pp-scattering changes its sign around ^/s ~ 30 GeV. At 
this energy the differential cross section of pp-scattering has its sharpest minimum. 

At the highest ISR energies all models and experiments show that p{0)pp > p{0)pp. The 
experiments also show that the dip at these energies in pp is higher then in pp-scattering. 

The observed dip-bump structure in the high-energy differential cross section was not 
predicted by any model or theory. It can be related e.g. to the multiple scattering of quarks 
and gluons (Glauber theory), leaving however much room for speculations on the basic inputs 
in this approach. The optical model (see Ref. [115] and earlier references therein) predicted 
a sequence of minima and maxima, however the ISR and SPS data show only one structure, 
both in pp and in pp scattering. Some models [116] predict that more structure will appear 
at higher energies, for example at the LHC. The (non) appearance of additional minima and 
maxima at the LHC will confirm or rule out part of the existing models, although most of 



them will survive by refitting its parameters a posteriori. 

Two models [22,117,118], fitted to the data in a wide range of s and t, can illustrate 
the state of art in this field (see, for instance. Fig. 8). While the authors of Ref. [22] claim 
that fits require the presence of an Odderon contribution, those of Ref. [117, 118] do not 
need it. By this we only want to stress that the flexibility of the existing models allow for 
good "postdictions" , but their predictions are not unique. In this situation useful empirical 
parameterization [119], unbiased by any theoretical prejudice can be useful. 
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Figure 8: Left panel: Elastic pp differential cross section predicted in Ref. [117,118]. Right 
panel: Elastic pp and pp differential cross section predicted in Ref. [22]. 

By definition, any diffractive pattern in da/dt is a property of the (predominantly imagi- 
nary) Pomeron contribution, rather than of the "non-diffractive" Odderon, whose amplitude 
is predominantly real and thus it can only "contaminate" the diffractive pattern. 

5.1 A simple model of the diffraction pattern 

Data on pp scattering below the LHC energy region can be described by a Pomeron and sub- 
leading contributions. From the phenomenological point of view, the need for the Odderon, 
becomes important only after the SPS data on pp scattering appeared, in which the dip seen 
in pp scattering appeared as a "shoulder" , suggesting that the diffractive minimum was filled 
by the Odderon contribution. This is not a proof, just a plausibility argument. Below we 
illustrate the dynamics of the dip and shoulder in pp and pp scattering by using a simple 
("minimal") dipole Pomeron (DP) model to "guide the eye". It reproduces the observed 
dynamics of the dip and will help to anticipate the LHC phenomena. 



By neglecting spin effects, the scattering amplitude for pp and pp scattering is written as 
a sum 

All{s,t)=Pis,t)±Ois,t), (23) 

where P is the Pomeron contribution and O is that of the Odderon. For the DP amplitude 
we use the following "minimal" model [23,24,26,27,120] 



P(s,t) = ^^(r2(s)e^?(^)["^W"i] - er2(s)e^'W["^W-^n, (24) 
bso ^ ^ 

where r^(s) = 6 + L — y, t\{^s) = L — ^ with L = In ap{t) is the Pomeron trajectory 
and a, 6, sq and e are free parameters. The model produces rising cross sections without 
violation of the Froissart bound as well as a dip-bump structure seen in the pp scattering 
with its non-trivial dynamics observed at the CERN ISR in the range 23 < -^i < Q2GeV 
[23,24,26,27,120]. The absence of a relevant structure (replaced by a "shoulder") in pp 
scattering at the CERN SPS was interpreted [23,24,26,27] as a manifestation of the Odderon, 
filling the dip produced by the Pomeron term (for reviews see Refs. [10,38]). A high-quality 
description of elastic high-energy data with predictions for future accelerators can be found 
in [116]. In an alternative successful approach [121] there are several Pomerons, whose 
interference produces the dip-bump structure and perfect agreement with the data. 

Since little is known about the properties of the Odderon, apart from its assumed asymp- 
totic nature, one usually parameterizes it in the form close to that of the Pomeron given in 
Eq. (24). However: 1) the Odderon is C— odd, which implies an extra i factor in front of 
the amplitude; 2) its relative contribution is by orders of magnitude smaller than that of the 
Pomeron (until now the Odderon was not seen in the forward direction); 3) the slope of the 
Odderon trajectory is much smaller than that of the Pomeron. There are two reasons for 
the latter: one is based on theoretical arguments [23, 24] and the other one is phenomeno- 
logical: the "flat" shoulder in pp scattering, seen at the SPS, could be a manifestation of the 
Odderon. 

Let us remind that the dip in the ISR energy region is monotonically deepening, reaching 
its maximal depth at 52.8 GeV, whereafter the monotonic trend changes, albeit at a single 
energy equal to 62 GeV. Similar to the unique case of the measured [34,35] difference between 
pp and pp amplitudes at 53 GeV, mentioned in the Introduction, this phenomenon needs 
confirmation. Regrettably, measurements of the difference between pp and pp amplitudes 
are not foreseen in the near future. 

Fig. 9 shows the differential cross section of pb scattering at 546 and 1800 GeV calculated 
from Eqs. (23) and (24) (full line), compared with a fit from paper [83] (broken line). Details 
of the calculations and the values of the fitted parameters can be found in Refs. [23,24,26, 
27,120]. 

It should be reminded that the hypothetical Odderon contributes enters Eq. (6) with 
different signs in pp and pp amplitudes, thus distorting the pure Pomeron contribution in 
both cases. Due to its small slope (flatness of its trajectory), the role of the Odderon increases 
with increasing — t, the phenomenon reaching its maximum in the dip region. 
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Figure 9: Differential cross section at 546 and 1800 GeV calculated from Eqs. (23) and (24) 
with the P + O contribution, corresponding to pp scattering (full line), compared with those 
of Ref. [83] (broken line). For detail of the calculations see Refs. [23,24,26,27,120]. 



The ratio of the cross section at the minimum to that at the maximum (depth of the 
dip) is more informative than its absolute value. Even greater than at the ISR (~ 5) value 
of the ratio will indicate the monotonic deepening of the dip and will disfavour the Odderon 
contribution, while its shallowing will favour the presence of the Odderon. 

At this point it may be appropriate to introduce several known models describing the 
diffractive pattern (the dynamics of the dip-bump structure). Before doing so, two points 
are worth mentioning: 1) one still lacks a theoretical understanding of the phenomenon and 
must rely on models; 2) among the large number of models, only few are able to fit properly 
the large number of high-statistics data that exist in a wide range of s and t for pp and 
pp scattering. A collection/review of the existing models for high-energy scattering with 
a critical evaluation of their prediction is highly desirable, however this task is beyond the 
scope of the present paper. We apologize to those authors whose papers, for brevity of space, 
are not cited here. 

5.2 The "Protvino" model 

One of the examples of eikonal models of diffraction is the "Protvino" model [92, 121]. 
The unitarity condition 

3mT(s,6) = \T{s,b)\^ + r]{s,b) , 



where T{s, h) is the scattering amphtude in the impact representation, h is the impact pa- 
rameter, ?7(s, h) is the contribution of inelastic channels, implies the following eikonal form 
for the scattering amplitude T(s, h): 

2i5{s,b) _ 1 

nsM = , (25) 

where h) is the eikonal function. The unitarity condition in terms of the eikonal looks 
as follows: 

3m 6) > 0, s > Sinei . (26) 

The eikonal function is assumed to have simple poles in the complex J-plane and the 
corresponding Regge trajectories are normally used in the linear approximation 

a(t) = a(0) + a'(0)t . (27) 

The following representation for the eikonal function is used: 

5%{s, h) = 5t,{s, h) + 5+ (s, h) + 5+ (s, h) T 5o(s, h) + 5j{s, b) TSZis,b), (28) 

where 6p_^^^{s,b) are Pomeron contributions. '+' denotes C even trajectories (the Pomeron 
trajectories have the following quantum numbers 0~^J~^~^), ' — ' denotes C odd trajectories, 
(5o(s, 6) is the Odderon contribution (the Odderon is the C odd partner of the Pomeron 
with quantum numbers 0~J ); 5^" , S~{s, b) are the contributions of secondary Reggeons, / 
(C = +1) and uj {C = -1). 

The parameters of secondary Reggeon trajectories are fixed according to the parameters 
obtained from a fit of the meson spectrum: 

af{t) = 0.69 + 0.84t , 

(29) 

a^{t) = 0.47 + 0.93t . 

The model fits high-energy elastic pp and pp scattering data. The data are well de- 
scribed for all momenta (0.01 < \t\ < 14. GeV'^) and energies (8. < ^/s < 1800 GeV) 
(xVd-o.f. = 2.60). It predicts the appearance of two dips in the differential cross-section 
which will be measured at LHC. The parameters of the Pomeron and Odderon trajectories 
are: 

a(0)p, = 1.058, a'{0)p, = 0.560 GeV^. 
a(0)p2 = 1.167, a'(0)p2 = 0.273 GeV-^; 
a(0)p3 = 1.203, a'(0)p3 = 0.094 GeV-^. 
a(0)o = 1.192, a'{0)o = 0.048 GeV-\ 

The model predicts (see Fig. 10) the following values of total and elastic cross sections 
at the LHC: 
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Figure 10: Left panel: Total cross sections of pp scattering (hollow circles) and pp scattering 
(full circles) and curves corresponding to their description in the three-Pomeron model [92, 
121]. Right panel: Predictions of the three-Pomeron model [92, 121] for the differential 
cross-section of pp scattering to be measured at LHC with -^i = 14. TeV and at RHIC 
y/s = 100. GeV and = 500. GeV. The data corresponding to the energy ^/s = 52.8 GeV 
is multiplied by 10-^ RHIC at 500 GeV by 10-^°, RHIC at 500 GeV by 10"^^ and that of 
LHC by 10-16. 



pp 

Otot 



106.73 (mb) 



+7.56 mb 
- 8.50mb 5 



(30) 



pp 

^ elastic 



29.19 (mb) 



+3.58 mb 
-2.83 mb ) 



ffP = 0.1378 



+0.0042 
-0.0061 • 



On the right panel of Fig. 10 one can see a typical "diffractive" pattern of the differential 
cross-section also present in other models based on unitarisation [117,118,122-125] (see also 
Fig. 9 and Fig. 11). 



5.3 The "Connecticut" model 



A model based on a physical picture of the nucleon having an external cloud, an inner shell 
of baryonic charge, and a central quark-bag containing the valence quarks was proposed 
in Ref. [93]. The underlying field theory model is the gauged Gell-Mann - Levy linear 
cr-model with spontaneous breakdown of chiral symmetry, with a Wess-Zumino-Witten 
(WZW) anomalous action. The model attributes the external nucleon cloud to a quark- 
antiquark condensed ground state analogous to the BCS ground state in superconductivity 
- an idea that was first proposed by Nambu and Jona-Lasinio. The WZW action implies that 




Figure 11: Left panel: Total cross sections of pp scattering (open circles) and pp scattering 
(full circles) and the fitted curves based on the model [93]. Right panel: Predictions of the 
model [93] for the differential cross-section of pp scattering to be measured at the LHC at 
^/s = 14. TeV (solid line) and comparison with the prediction of other models: Avila et 
al. [22], Block et al. [126], Bourrely et al. [117,118], Desgrolard et al. [127], and Petrov et al. 
(three pomerons) [92]. 



the baryonic charge is geometrical or topological in nature, which is the basis of the Skyrmion 
model. The action further shows that the vector meson u couples to this topological charge 
like a gauge boson, i.e. like an elementary vector meson. As a consequence, one nucleon 
probes the baryonic charge of the other one via cj-exchange. In pp elastic scattering, in the 
small momentum transfer region, the outer cloud of one nucleon interacts with that of the 
other giving rise to diffraction scattering. As the momentum transfer increases, one nucleon 
probes the other one at intermediate distances and the a;-exchange becomes dominant. At 
momentun transfers even larger, one nucleon scatters off the other via valence quark-quark 
scattering. 

Diffraction is described by using the impact parameter representation and a phenomeno- 
logical profile function: 

Tn{s,t)=ipW bdb Joibq)TDis,b) . (31) 
Jo 

here q is the momentum transfer (g = ^J\t\) and TD{s,b) is the diffraction profile function, 
which is related to the eikonal function xo(s, 6): T£){s,b) = 1 — exp{ixD{s,b)). T£){s,b) is 
taken to be an even Fermi profile function: 



^D{s,b) = g{s) 



+ 



1 + (,{b-R)/a I ^-{h+R)/a 



- 1 



(32) 



The parameters R and a are energy dependent: 

R = Ro + Ri (\n—-i- 
V So 2 

a = ao + fli ( In — - ) ; 

V So 2/ 

g{s) is a complex crossing even energy- dependent coupling strength, so = 1 GeV^. 
For the diffraction amplitude the following asymptotic properties are obtained: 

1- o"tot(s) ~ (flo + flilns)^ (Froissart-Martin bound) 

2. p(s) ^ ao+ai\ns (derivative dispersion relation) 

3. T£){s,t) i s ln^s/(|t| In^ s) (Auberson-Kinoshita-Martin scaling) 

4. T^{s,t) = T^{s,t) (crossing even) 

At 14 TeV atot and ppp were found to have the values 110 mb and 0.120 mb, respectively. 

Finally, Fig. 11 shows a compilation of predictions of the model for pp elastic dcr/dt at 
LHC together with the predictions of other models, e.g.: Avila et al. [22], Block et al. [126], 
Bourrely et al. [117,118], Desgrolard et al. [127], and Petrov et al. (three pomeron) [92]. 
As one can see from Fig. 11, the model [93] predicts no diffractive pattern at large \t\ as 
diffractive scattering dominates in the small \t\ region, hard scattering dominates in the 
intermediate |t| region (1.5 < \t\ < 8.0 GeV^) and quark-quark scattering takes over at large 
\t\. 



5.4 The "Dubna Dynamical" model 

The so-called "Dubna Dynamical" (DD) model was already cited in connection with the 
nearly forward scattering; it will be used also in Sec. 6.2 on the BDL, so we find it appropriate 
to briefly introduce it, especially as it may be closely related to models proposed earlier (see 
Refs. [117,118,122-124]) or later (see Refs. [92,93,121] and earlier references therein). 

The DD model was developed in Refs. [125,128,129] to describe particle interaction with 
account for the internal structure of hadrons. According to the DD model, the hadron 
scattering amplitude can be written as a sum of a central and a peripheral part of the 
interaction [125]. The full scattering amplitude is calculated using the impact parameter 
representation with the eikonal form of the unitarization: 

A{s,t) = is / bdbJo(bA)[l - 6-"^^''% (33) 
Jo 



where 



X{s, b) = Xc(s, b) + Xp(s, b) + xr{s, b). 



(34) 



Here the central interaction term Xc{s, b) describes the interaction between central parts of 
the hadrons. At high energies, it is determined by a spinless Pomeron exchange. The second 
term, Xpi^i b) is the sum of triangle diagrams corresponding to the interaction of the central 
part of a hadron with the meson cloud of the other one. The meson- nucleon interaction leads 
to spin flip effects in the Pomeron-hadron vertex [114]. The term Xr{s, b) describe the cross 
even and cross odd parts of the subleading Reggeons with the energy dependence ~ 1/ ^/s. 



Figure 12: Predictions of the DD model [88,89] for the Tevatron energies, ^/s = 1.8 TeV 
(long dashed line) and for the LHC energies (short dashed line). For the sake of comparison, 
model curves and experimental data are shown for one lower energy, ^/s = 23 GeV (solid 
line). The experimental data points are from Ref. [130]. 

The term Xc{s, b ) was chosen phenomenologically (for details see Refs. [122-124]), where 
the scattering amplitude was obtained from an integral representation of the Bessel function. 
It was shown in those papers that, starting from the usual partial wave expansion in the 
Legendre polynomials, the eikonal representation for the scattering amplitude can be derived 
to be valid for all energies and in the whole angular domain, and for which the spectral 
function is given in terms of the partial wave amplitudes. The scattering amplitude can be 
presented in the form 



Although the analytic calculations can be performed to any u, to simplify numerical calcu- 
lations, we set = 1/2, whereafter the amplitude becomes: 





(35) 




(36) 



By taking this form as the Born amphtude, the eikonal term xo in the DD model assumes 
the form 

Xo(s,6) = e-™W (37) 

It leads to an exponential behavior of the eikonal at large values of the impact parameter 
and a Gaussian behavior at small impact parameters. In the DD model the inelastic states 
in the s-channel were also taken into account. Accordingly, the central part can be presented 
as 

Xci-s, b) = xoi-s, b) [1 - 7Xo(s, b)], (38) 

with 

Xois,b) = Ao e-""^'^ (39) 

The coefficient 7 is connected with the contribution of the inelastic intermediate state in the 
s-channel. Here Ao, m(s) and Rq are free parameters, whose values and energy dependence 
were determined from fits to the experimental data on proton-proton elastic differential 
cross section scattering from ^/s = 9.78 GeV to ^/s = 540 GeV (see Fig. 12) and momentum 
transfer from t = to t = —12 GeV^. 

The additional term Xp determined by the peripheral Pomeron interaction with 
account for the meson cloud. This leads to the important feature of the DD model, namely 
the existence in the eikonal phase of a small term growing as ^/s : 



Xp{s, b) = H,o + Av^ rrip^Rl + b^K^{ m,^ Rl + h^), (40) 

where Ki is the MacDonald function of the first order. So, the contribution to the eikonal 
phase growing as -y/i is determined by the peripheral meson-cloud effects. It has been shown 
that this term becomes important for energies \/s > 200 GeV. The corresponding form of 
the interaction potential can be represented, see Ref. [122], by a superposition of Yukawa- 
like potentials in the form of the MacDonald function. The leading term of the Xc in the 
momentum transfer representation is 



~ (41) 



/4m2 + g2 

It results in the following form of the differential cross section for » [125] 



dt t 



2 



(42) 



Notice that the square-root behavior in the above formula corresponds to the square-root 
Pomeron trajectory, Eq. (2), both following from analyticity and unitarity, with important 
physical consequences, such as the "break" in the cone near t = —0.1 GeV^ [10,39] or levelling 
off of the impact parameter amplitude at large b, corresponding to the mesonic atmosphere 
of the nucleon. 

The model provides a self-consistent picture of the differential cross section and spin 
phenomena of different hadron processes at high energies. Really, the parameters in the 
amplitude determined from one reaction, for example, elastic pp-scattering, allow to give 
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Figure 13: Predictions of the DD (solid and dashed hnes) and BSW models (triangles and 
full squares) for the differential cross section, at y/s = 2 TeV (dashed line and triangles) and 
a/s = 10 TeV (solid line and full squares). 



predictions on elastic meson-nucleon scattering and charge-exchange reaction 7r"p n^n 
at high energies. The model predicts that at super-high energies polarization effects for 
particles and antiparticles are the same [114,131]. 

Let us note that a similar model, describing simultaneously different hadronic elastic 
cross sections, polarization and spin correlation effects was developed by Bourrely, Soffer 
and Wu [117, 118] (BSW). According to the physical picture of this model, each hadron 
appears as a black disk with a gray fringe, where the black disk radius increases as ln(s). It 
is interesting to compare the predictions of the BSW and DD models for the LHC energies. 
The result is shown in Fig. 13 and Table 6. The most important difference comes from the 
region \t\ > 0.3 GeV^. 

Future experiments at the LHC will provide for an excellent possibility to test the theo- 



energy 


model 


CTtot, mb 


CTeh mb 


crtnt 




Pit = 0) 


f mb/GeV^ 


2 TeV 


DD 

BSW 


81 
76.1 


20.7 
17.9 


0.256 
0.235 


83 
82 


0.197 
0.128 


38.8 
11.2 


10 TeV 


DD 
BSW 


123 
98.4 


42.6 
26.8 


0.35 
0.27 




0.195 
0.122 


350 
38.6 



Table 6: Comparison of the predictions from the DD and BSW models for two LHC energies. 



retical argument on the peripheral rise of the eikonal phase at superhigh energies determined, 
in this approach, by the meson-cloud effects. 



6 Intermediate and Large \t\ 

The definition of the transition between "soft" and "hard" physics is rather ambiguous, 
although there are phenomena indicative of this transition. One is the slow-down of the 
shrinkage of the cone, already noticed at the ISR, and attributed to the hard scattering of 
the constituents of nucleons. It will manifest also in the change of the exponential decrease to 
a power fall-off in Furthermore, in this region the exponential decrease of the scattering 
amplitude will be gradually replaced by a power one. These effects have a simple interpreta- 
tion within a Regge-pole model with nonlinear trajectories. As shown in a number of papers 
(see e.g. Refs. [57,59]), logarithmic trajectories in the dual Regge-model [55,60] can mimic 
hard scattering of point-like constituents and the transition from the soft interaction of the 
extended objects, like strings, to the hard scattering of point-like particles, governed by the 
quark counting rules and/or QCD. 

alpha' 




t 



Figure 14: Typical behavior of the slope of a logarithmic trajectory. Using Eq. (44), the 
slope a'{t) decreases between t = —1 and t = —10 by about a factor 3. According to Eq. (43) 
the shrinkage of the slope is expected to slow down by the same rate. 

A particularly simple expression for the slope of the diffraction cone, illustrating the 
aforesaid phenomena, follows from the geometrical properties of the geometrical model (see 
Eq. 24 and Ref. [132]): 

B{s,t) = ka'{t)at{s), (43) 

where /c is a parameter and a'{t) is the slope of the logarithmic Pomeron trajectory given 
by Eq. (3). By choosing 7 = 1 and P = 0.3, we get from Eq. (3) 



// N 0.3 

a'(t) = . 

^ ^ 1 - 0.3t 



(44) 



The slope a'{t) given by Eq. (44) is shown in Fig. 14. The logarithmic trajectory (3) 
(or, similarly, (4)) mimics also the transition from the exponential decrease of the scattering 
amplitude in t (soft physics) to a power one, according to 



[b+\n{s)] In(-t) 



(45) 
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Figure 15: Predictions of the DD model for the differential cross section at ^/s 
(solid line) and at = 14 TeV (dashed line) [88, 89] 



52 GeV 



Note that in this region the ratio \t\/s is much smaller (tends to zero) than at the ISR 
energies. The new information obtained in this region (small angles and large t) will give the 
possibility to determine whether perturbative QCD works here. The logarithmic regime of 
the trajectory in the above Regge-pole model with the resulting power-like behavior in this 
kinematical region provides a link with the quark-parton model and, eventually, with the 
perturbative regime of QCD. However this transition is much more complicated, involves a 
rearrangement of the dynamics, that can be treated e.g. in the dual model of Ref. [59]. 



The valence qq scattering is assumed to be due to a color singlet amplitude behaving as 
s/t, so that daqq/dt ~ For large this leads to a power-like behavior [133]: da/dt ~ 
t"^*^, which corresponds to the dimensional counting rules [134,135]. The perturbative QCD 
behavior can be preceded by a transition regime da/dt ~ t~^. 

Different models of "hard" scattering, with power-like scattering amplitudes obeying 
scaling, exist in the literature. For example, in a quark-diquark model of the nucleon [136], 
the introduction of the transverse distance 1 / ^J\t\ between the colliding valence quarks gives 
an additional degree of freedom. 

Hard scattering and the power-low behavior of the cross sections by itself is an important 
issue of the strong interaction dynamics, which, however goes beyond the scope of the present 
paper. In large-t elastic pp scattering at the LHC one will enter only the transition region 
between the "soft" (exponential) and "hard" power-like behavior of the cross section. 

Contrary to Eq. (45), that predicts the decrease of the differential cross section in this 
region, some models lead to an increasing elastic differential cross section at large momentum 
transfer. For example, the DD model [88, 89] describes the elastic differential cross section 
not only at small \t\, but also up to large \t\ ~ 10 GeV^ at the ISR (see Fig. 15). Note that 
the experimental data show that, up to the TeV-region, the differential cross section has a 
weak energy dependence at large t. This effect should be better understood, since in the 
TeV-region a change of the regime might occur, as Fig. 15 seems to indicate. 



7 Black disc limit at the LHC? 



There is general consensus on the hypothesis that with increasing energy nucleons are becom- 
ing blacker, edgier and larger (BEL effect). Since the details depend on model calculations 
and their extrapolations beyond the energies of present accelerators, the onset and conse- 
quences of BEL effect also differ. Of interest is the possibility that at the LHC the nucleons 
will reach the black disc limit (BDL) with measurable effects. If or when the black limit will 
be reached (at the LHC?) two cases, connected with the use of the s— channel unitarity, are 
possible. According to the first option, the BDL is identical to the unitarity limit, absolute 
for the central opacity of the nucleon. According to the second option (an alternative solution 
of the unitarity equation, which will be illustrated below), the unitarity limit is well beyond 
the BDL, leaving space for a further increase with energy of the central "overlap function" 
(or impact parameter amplitude). In this approach, the nucleon, after having reached the 
BDL and tending to the unitarity limit, will gradually become more transparent. Below we 
present both options and their observable consequences. 



7.1 Definitions 

The onset of "saturation", wherever it occurs, can hardly be seen directly from the data, 
since the behaviour of the real part of the amplitude is at best known in the nearly forward 
direction only, while for the rest it can be only modelled. The expected approach to the BDL 



depends both on the model for the scattering amphtude which is used and on the procedure 
of unit arizat ion. Numerical estimates based on particular model calculations are presented 
in the two subsequent subsections. 

To start with, let us remind the general definitions and notations. Unitarity in the impact 
parameter b representation reads 

m{s,b) = \h{s,b)f + G{s,b). (46) 

Here h{s, b) is the elastic scattering amplitude at the center of mass energy ^/s, Qh{s, b) is 
usually called the "profile function", representing the hadron opacity and G{s,b), called the 
"inelastic overlap function", is the sum over all inelastic channel contributions. Integrated 
over 6, Eq. (46) reduces to a simple relation between total, elastic and inelastic cross sections: 

Eq. (46) imposes the absolute limit 

<\h{s,b)f <'^h{s,b) <1, (47) 
while the so-called black disc limit 

aeiis) = ainis) = ^atotis) 

or 

Qh{s, 6) = 1 (48) 

is a particular realization of the optical model, namely it corresponds to the maximal ab- 
sorption within the eikonal unitarization, when the scattering amplitude is approximated 
as 

his,b) = '-il-exp[iiuis,h)]), (49) 
with a purely imaginary eikonal ci;(s, 6). 

The solutions of the unitarity equation (46), with ^h{s,b) = 0, are 

1 



h{s,b) 2 



1 ± Jl-AG{s,b) 



(50) 



The solution with the minus sign in front of the square root describes the eikonal unitariza- 
tion. The other one, that with the plus sign in front of the square root, is known within the 
so-called ?7-matrix approach [132,137-139], where the unitarized amplitude has the form of 
a ratio rather than of an exponential, typical of the eikonal approach. It is written as 

where U is the input Born term, the analogue of the eikonal u in Eq. (49). The connection 
between the two forms was discussed in recent papers [140,141]. 



In the [/-matrix approach, the scattering amplitude h{s,b) may exceed the black disc 
limit as the energy increases. The transition from a (central) black disc to a (peripheral) 



black ring, surrounding a gray disc, for the inelastic overlap function in the impact parameter 
space corresponds to the transition from shadowing to antishadowing [139] . We shall present 
a particular realization of this regime. 

The impact parameter amplitude h{s, b) can be calculated either directly from the data, 
as it was done e.g. in Ref. [142] (where, however, the real part of the amplitude was ne- 
glected) or by using a particular model that fits the data sufficiently well. We consider three 
representative examples, namely the Donnachie-Landshoff (DL) model [1,2,143-146], the 
DP model [26,27,132] and the DD model [114,125,128,129]. 



7.2 The "Born term" 

The Donnachie-Landshoff model [1,2,143-146] is very popular for its simplicity. Essentially, 
it means the following four-parametric empirical formula to fit all total hadronic cross sections 

atot = X s' + Y s'^, (52) 

where two of the parameters, namely 6 = ap{0) — 1 ^ 0.08, and 5r, which is negative, are 
universal. The violation of the Froissart-Martin (FM) bound, 

atotis) <C {In sY, C = 60mb, (53) 

inherent to this model, is rather an aesthetic than a practical defect. 

The t dependence in the DL model is usually chosen [143-146] in the form close to the 
dipole Pomeron form factor. For the present purposes a simple exponential residue in the 
Pomeron amplitude will do as well, with the signature included: 



A(s,t) = - N i-i—] e^* , (54) 
V Sdl J 

where a(t) = ap(0) +a' t is the Pomeron trajectory and is a dimensionless normalization 
factor related to the total cross section at s = Sdi by the optical theorem: 

^ = ~A ~^ 7m ^i"**^^ = • ^^^"f 

47rsm |ap(0) 

According to fits of Ref. [143-146] one has sai = 1 GeV^, ap(0) = 1.08, a' = 0.25 GeV-^, 
and X = 21.70 mb (see (Eq. 52)) and thus N = AAA. By puting 

Ms,t) ^ da{s,t = 0) 

dt dt ^ ' 

and choosing the CDF or E410 results for the slope -Bexp at the Tevatron energy, we obtain 
B = 5e,.p(s)/2 - a'\n{s/sdi) = 4.75 GeV^. 



In the DP model [132], factorizable at asymptotically high energies, logarithmically rising 
cross sections are produced at the Pomeron intercept equal to one, so that the DP model 
does not confiict with the FM bound. While data on total cross section are compatible with 
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a'(GeV-2) 


e 


so(GeV2) 


355.6 


10.76 


1.0356 


0.377 


0.0109 


100.0 



Table 7: Values of the parameters of the DP model, Subsections 5.1 and 7.2, found in 
Ref. [26,27]. 

a logarithmic rise (according to the DP model), the ratio (yf,ijotot is found (see Ref. [132] for 
details) for 5 = to be a monotonically decreasing function of the energy for any physical 
value of the parameters. The experimentally observed rise of this ratio can be achieved only 
for 5 > and thus requires the introduction of a supercritical Pomeron, with a(0) > 1. As 
a result, the rise of the total cross section is driven and shared by the dipole Pomeron and 
the "supercritical" intercept. The parameter 5 = ap(0) — 1 in the supercritical DP model 
is nearly half that of the DL model, making it safer from the point of view of the unitarity 
bounds. Generally speaking, the closer the input to the unitarized output, the better the 
convergence of the unitarization procedure. 

Let us remind that, apart from the conservative FM bound, any model should satisfy 
also s-channel unitarity. We demonstrate below that both the DL and DP models are well 
below this limit and will remain so for long time, in particular at LHC. Let us remind that 
the DL and the DP model are close numerically, although they are different conceptually 
and consequently their extrapolations to superhigh energies will differ as well. 

The elastic scattering amplitude corresponding to the exchange of a dipole Pomeron reads 

= e-™/2(^)"[G'(a) + (L-z7r/2)G(a)], ^ ^ 



where L =ln^ and a = ap{f) is the Pomeron trajectory. 



By setting G'[ol) = —ae'^^^" ^\ Eq. (57) can be rewritten in the geometrical form dis- 
cussed in Subsection 5.1, (see Eq. (24)). 

In Table 7 we quote the numerical values of the parameters of the DP model fitted in 
Ref. [26, 27] to the data on proton-proton and proton-antiproton elastic scattering: 

crtotis) = ^QA{s, 0) , p{s) = 1^1^ ; 4 (GeV) < < 1800 (GeV) (58) 

as well as the differential cross-section 

^^ipH = ^M(s,t)|',23.5(GeV)<v^< 630 (GeV),0(GeV')<|t|<6(GeV2). (59) 
dt 

In that fit, apart from the Pomeron, the Odderon and two subleading trajectories u and / 
were also included. Here, for simplicity and clarity we consider only the dominant term at 
high energy due to the Pomeron exchange with the parameters fitted in Ref. [26,27]. 



We use the above set of parameters to calculate the impact parameter amplitude. 





53 GeV 


546 GeV 


1800 GeV 


exp 
th 


0.36 
0.36 


0.420 ± 0.004 
0.424 


0.492 ± 0.008 
0.461 



Table 8: Central opacity of the nucleoli '^h{s, 0) calculated at ISR, SPS and Tevatron energies 
compared with experimental results. 

7.3 Impact parameter representation, the black disc limit 
and unitarity 

The elastic amplitude in the impact parameter representation in our normalization is 

1 

h{s,b) = — dq qJo{bq)A{s,-q^), q = , (60) 

2s Jo 

expressed in tems of the Bessel function Jq. 

The impact parameter representation for linear trajectories is calculable explicitly for 
the DP model using Eq. (57). We have 

h{s,b) = I go [e'-"' e-'""^"^ - e e'"' e-^'l^^l\ , (61) 
where the functions rl = rf{s){i = 1, 2) have been introduced in the Subsection 5.1 and 

i?? = «V^ 9o = jr^. (62) 

Asymptotically (i.e. when L ^ bp, which implies ^/s ^ 2. TeV, with the parameters of 
Table 7) we get 

h{s, b) i g{s) {1 — e) e ^ , s — > oo, (63) 
where ^ 

R^ = a'L, g[s) = go{^j) . (64) 



To illustrate the effect coming from s-channel unitarity, in Ref. [148, 149] a family of 
curves showing the imaginary part of the amplitude in the impact parameter-representation 
as well as the inelastic overlap function G(s, b) calculated from Eq. (46) at various energies, 
was displayed. 

Our confidence in the extrapolation of Qh{s, b) to the highest energies rests partly on 
the good agreement of our (non fitted) results with the experimental analysis of the central 
opacity of the nucleon (see Table 8). 

It is important to note that the unitarity bound 1 for Qh{s, b) will not be reached at the 
LHC energy, while the black disc limit 1/2 will be slightly exceeded, the central opacity of 
the nucleon being Qmh{s, 0) = 0.54. 



^Similar calculations with non-linear trajectories can be found e.g. in Rcfs. [132,147]. 



The black disc limit is reached at ^/s ~ 2 TeV, where the overlap function reaches its 
maximum value 1/4. This energy corresponds to the appearance of the antishadow mode in 
agreement with the general considerations in Ref. [139]. Notice that while '^h{s,b) remains 
central all the way, G{s, b) is getting more peripheral as the energy increases starting from the 
the value of Tevatron. For example at ^/s = 14 TeV, the central region of the antishadowing 
mode, obtained from the U matrix unitarization, below b ~ 0.4 fm is discernible from 
the peripheral region of shadowing scattering beyond b ~ 0.4 fm, where G{s,b) = 1/4. 
Consistently with the BEL effect the proton will tend to become more transparent at the 
center (in the sense of becoming a gray object surrounded by a black ring), i.e. it is expected 
to become gray, edgier and larger (GEL effect). 

The s channel unitarity limit will not be endangered until extremely high energies (10^ 
Gev for the DL model and 10^ GeV for the DP model), safe for any credible experiment. It 
is interesting to compare these limits with the one imposed by the FM bound: actually the 
Pomeron amplitude saturates the FM bound at 10^^ GeV. As expected, the FM bound is 
even more conservative than that following from s-channel unitarity. 

Now, we consider the unitarized amplitude according to the tZ-matrix prescription [132, 
137, 138] 

H{s,b)= , (65) 

1 — th[s, b) 

with the Born term h{s,b) calculated from Eq. (60) using Eq. (57). 

An unescapable consequence of the unitarization is that, when calculating the observ- 
ables, one should also replace the Born amplitude A{s, t) with a unitarized amplitude A{s^ t) 
defined as the inverse Fourier-Bessel transform of h{s, b): 

A(s,t) = 2s dbb.Ub^)His,b) . (66) 

JO 

Thus, the above picture may change since the parameters of the model should in principle 
be refitted under the unitarization procedure (this effect of changing the parameters was 
clearly demonstrated e.g. in Ref. [150]). 

While the unitarity limit now is secured automatically (remind that '^h{s, 0) is well below 
that limit even at the Born level in the TeV region), the behaviour of the elastic impact pa- 
rameter amplitude after it has reached the black disc limit corresponds (see Ref. [139]) to the 
transition from shadowing to antishadowing. In other words, the proton (antiproton) after 
having reached its maximal blackness around 2 TeV, will become gradually more transparent 
with increasing energies at its center. It follows from the presented model that, in getting 
edgier and larger, the proton, after reaching its maximal blackness, will tend to be more 
transparent (GEL effect), i.e. a gray disc surrounded by a black ring will gradually develop 
beyond the Tevatron energy range. The transition from shadowing to a new, antishadowing 
scattering mode is expected to occur at the LHC. 

To conclude, we stress once more that both the data and relevant models at present 
energies are well below the s-channel unitarity limit. Deviations due to the diversity of 
realistic models may result in discrepancies concerning '^h{s, 0) of at most 10%, while its 
value at 2 TeV is still half that of the unitarity limit, so there is no reason to worry about 
it. Opposite statements may result from the confusion with normalization. Therefore, the 



p (v/i = 14 TeV, t) 



t = GeV 2 t = -0.1 GeV ^ 

DD [88,89] Soft and hard Pomeron [90] DD Soft and hard Pomeron 
0.19 0.24 0.08 0.05 



Table 9: The effect of saturation on the ratio p{t) as predicted in Refs. [88,89], [90]. 

model amplitudes at the Born level may still be quite interesting and efficient in analyzing 
the data at energies of present accelerator and giving some predictions beyond. The question, 
which model is closer to reality and meets better the requirements of the fundamental theory 
remains of course topical. 



7.4 Saturation at the LHC (in the DD model) 

The effect of the Black Disc Limit (BDL) saturation can be observed in the change of the 
i-dependence of the slope B and of the ratio p beginning from approximately ^/s = 2^6 TeV 
and tending to enhance at y/s = 14 TeV [90]. Such a saturation effect was studied in the 
framework of the DD model. 

In Refs. [90, 151] the BDL is determined by account of a hard Pomeron. The model 
predicts the increase with \t\ of the slope B{t) at small t and energies typical of the LHC, 
as shown in Fig. 16, right panel. In the same Figure, left panel, it is shown that saturation 
of the BDL substantially modifies the t-dependence of the ratio p{s,t) (see, also. Table 9). 
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Figure 16: Left panel: The ratio of the real to the imaginary part of the amplitude as a 
function of t, for the bare and the saturated amplitudes at various energies: 100 GeV (solid 
line), 500 GeV (long dashes), 5 TeV (short dashes) and 14 TeV (dash-dotted curve). Right 
panel: The slope of the elastic differential cross section as a function of t, for the bare and 
saturated amplitudes at various energies: 100 GeV (solid line), 500 GeV (long dashes), 5 
TeV (short dashes) and 14 TeV (dash-dotted curve). 



In fact, saturation naturally predicts a small increase of the slope B{t) with t at small t. 



To understand this phenomenon, let us take a simple black disk form with a sharp edge at 
radius R. The relevant scattering amplitude can then be written as 



A{s,t = 0) 



In this model the slope of the differential cross section at small momentum transfer 
behaves as 



Hence the slope grows with increasing |t| at small momentum transfer, as one can see in 
Fig. 16. It is interesting that two quite different models lead practically to the same result. 

Saturation is usually attibuted to the growing gluon density at small x. This phenome- 
onon is important at relatively large impact parameters, where the role of non-perturbative 
effects [152], due to confinement, is not clear. On the other hand, the BDL-saturation is 
generally connected with the small impact parameter region. Scrutinizing the small \t\ re- 
gion and looking for possible saturation effects there [153] is of great interest for future LHC 
measurements. 



8 Summary and outlooks 

In the present paper we concentrated mainly on elastic scattering and total cross sections. 
It should be remembered that the elastic scattering amplitude is the central, in a sense, 
basic object of the theory, therefore its knowledge is essential for the further progress in 
understanding strong interactions. Furthermore, elastic scattering, by unitarity (see, e.g. 
Sec. 7), is directly connected to inelastic processes. 

The LHC facilities provide a unique experimental environment where one can test prop- 
erties of the Pomeranchuk trajectory. Let us remind the open questions and problems that 
can be addressed both theoretically and experimentally: 

• The Pomeron intercept. From perturbative QCD the value of about a(0) ~ 1.25 follows 
[17-20]. It should be remembered, as often stressed by the authors of the BFKL equation, 
that this value should not be confronted with the data prior the unitarization of the amplitude 
(that will lower the rate of increase of the total cross section). A complete unitarization 
procedure, however, is difficult to be realized. In any case, experimental measurement of the 
total cross section will shed new light on the value of the Pomeron intercept. 

• Similar arguments concern the slope of the trajectory, which, because of its non-linear 
nature, is t-dependent. Hence it is important to explore the small \t\ region of the elastic 
differential cross section. 

• By parameterizing the diffraction cone by the method of overlapping bins, the fine 
structure of the cone (or that of the Pomeron) will be detected (Sec. 3). 




(67) 



• Perturbative QCD predicts logarithmic asymptotics for the Pomeron trajectory, as in 
Eq. (3). This basic result of QCD is often ignored. At the LHC, the onset of, or approach 
to the logarithmic behavior of the Pomeron trajectory should be visible. 

• The "QCD Pomeron" [17-20] implies an infinite number of poles or a moving cut, 
which are difficult to be realized phenomenologically. In a simplified version, they can be 
approximated as a sum of two poles with different intercepts. 

• The shape and the depth of the expected diffractive minimum and its ratio to the height 
of the second maximum are extremely informative (Sec. 5). The theoretical predictions here 
are not unique, but future data from LHC will narrow the margin for existing theoretical 
models. 

• The expected transition from soft to hard physics will be revealed at the LHC. The 
determination of the s dependence of the elastic cross section for |t| > 1 GeV^ will be very 
important for understanding the onset of the new, "hard" dynamics, replacing the "soft" , 
diffractive one. 

• A basic characteristics of the strong interaction dynamics is the ratio R{s) = cjei/uti 
reflecting the intensity of particle production and/or the realization of a particular unitariza- 
tion scheme. In the past, R{s) was changing unpredictably: from the (seemingly asymptotic) 
constant in the ISR energy region, suggesting geometrical scaling, to a subsequent indeflnite 
rise. Will this trend continue or a "new asymptotic" plateau is expected at the LHC? 
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